FLOWFIELD ANALYSIS FOR SUCCESSIVE OBLIQUE SHOCK WAVE INTERACTIONS 
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. A computation procedure is described for predicting flow fields which 
develop when successive interactions between oblique shock waves and a 
turbulent boundary layer occur. Such interactions may occur, for example, 
in engine inlets for supersonic aircraft. Computations have been carried 
out for axis ymme trie internal flows at = 3.82 and 2.82. The effect of 

boundary layer bleed has been considered for the = 2.82 flow. A control 
volume analysis is used to predict changes in the flow field across the 
interactions. Two bleed flow models have been considered. A turbulent 

i 

boundary layer program has been used to compute changes in the boundary layer 
between the interactions. The results given are for flows with t\i7o shock wave 
Interactions and for bleed at the second interaction site. In principle 
the method described may be extended to account for additional interactions. 
The predicted results are compared with measured results and are shown to be 
in good agreement when the bleed flow rate is low (on the order of 3 percent 
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of the boundary layer mass flow) , or when there is no bleed. As the 
bleed flow rate Is increased, differences between the predicted and 
measured results become larger. Shortcomings of the bleed £low models 
at higher bleed flow rates are discussed. 
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a constant in the wall-wake profile 

{a+{(Y-l)/2]M^)/(T„/V> - 1 

a constant In the Law of the Wall (usually equals 5.1) 

2 

skin friction coefficient, t^/U/2) 

X“in.onientiini of the bleed flow 

a constant in the Law of the Wall (usually equals O.A) 
shock wave-boundary layer interaction length 
Mach nuniber 

boundary layer mass bleed rate 
pressure 

radial coordinate from tunnel centerline 

radial coordinate of dividing stream surface separating 
bleed flow from main flow, see Figure 3. 

Reynolds Nuniber 

velocity in streamwise direction 
van Driest’ s generalized velocity, 

(u /A) arc sin {[(2A^ u/u^) - B]/(B + AA ) } 

^ ‘ft 

friction velocity, 

axial coordinate, measured from shock generator tip 
coordinate normal to the tunnel wall 
ratio of specific heats 

a thickness of freestream flow to allow for boundary layer 
mass entrainment, see Figure 3 

' boundary layer thickness 
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A 

6 displacement thickness of the boundary layer 

n “ y/5 


e 

V 

n 

p 


- momentum thickness of the boundary layer 
e kinematic viscosity 
= coefficient of the wake function 
*= mass density 


o “ I(Y-l)/2] M^/ <1 + t( -D/2] M^} 

Subscripts 

r> 

e conditions at the edge of the boundary layer 
w conditions at the wall 

*® frees tream conditions ahead of the first interaction 
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INTRODUCTION 

The Interaction of an oblique shock wave with a turbulent boundary 
layer is known to induce drastic changes in the boundary layer properties 
and to cause substantial deviation of the supersonic flow field from the 
predicted inviscid flow. This deviation may be of sufficient magnitude to 
adversely affect the performance of aerodynamic devices. Suitable methods 
for predicting the boundary layer and frees tream flow characteristics in 
the presence of such disturbances are required by engineers responsible 
for the design of aerodynamic configurations involving shock wave boundary 
layer interactions. 

A control volume method developed by Seebaugh, Paynter and Childs^, 

2 

and improved upon by Mathews » has been successful in the prediction of 
boundary layer characteristics downstream of the interaction with a single 
oblique shock wave- However, in some aerodynamic devices, such as mixed 
compression supersonic diffusers, the turbulent boundary layer is often 
subjected to interactions with more than one shock wave. In the study 
reported here the control volume method has been used^ to calculate the changes 
in turbulent boundary layer characteristics across successive shock wave 
interactions at the walls of axlsymmetric wind tunnels. Two cases are 
reported, one at « 2 . 82 , the other at M^ = 3 . 82 , In the Mach 2.82 study 
the effect of boundary layer bleed at the second interaction site was 
considered. A turbulent boundary layer computer program was used to predict 
boundary layer changes bet\-7een the two interactions and to provide the initial 
conditions for the second interaction. 

EXPERIMENTAL CONFIGURATIONS 

The experimental configurations which were used to produce the successive 
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shock waves are shown schematically in Figure 1. The Hach 3.82 tunnel 
had a radius of 1.015 inches and a boundary layer thickness ahead of the 
first interaction of 0.170 inches. The shock wave generator was installed 
on the centerline of the tunnel at zero angle of attack. The generator 
had a 10-degree half-angle conical tip which broke to 13 degrees 2.60 inches 
downstream of the tip. For the Mach 2.82 tunnel the tunnel radius was 1.02 
inches, and the boundary layer tliickness ahead of the first interactxon was 
0.165 inches. -The shock generator had a 10-degree conical tip which broke 
to 13 degrees 1.942 inches behind the tip. Both generators were designed 
to provide as large a region of freestream flow as possible between the 
first reflected and second Incident shock waves while at the same time keeping 
the expansion wave off the downstream comer of the second conical surface 
from interfering with the second interaction. 

^out 3.0 inches from the tip of the Mach 2.82 shock generator, or 
approximately at the pdnt where the second incident shock wave reached the 
wall, two rows of thirty-eight 0.052-inch diameter bleed holes were drilled 
around the periphery of the tunnel, ^ The bleed system‘was operated in a 
choked flow condition. With one row of the holes open, the bleed mass flux 
was about 2.8 percent of the boundary layer mass flux just ahead of the 
second interaction. With two rows of holes open, the bleed mass flux was 
5.0 percent of the boundary layer mass flux. 

Both tunnels were operated with a steady supply of dry air at 540 ®R. 

The freestream unit Reynolds number for the M *= 3.82 tunnel was 5,62 x 10 
per foot, that for the M “ 2.82 tunnel, 5.8 x 10 per foot* 
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INSTRUMENTATION 

Standard instrumentation was used to obtain tunnel wall static 
pressures and boundary layer pitot pressure profiles. Wall static pressures 
were taken at 0.05-incli intervals along the tunnel sidewalls. Pitot profiles 

r 

were taken in radial increments of 0.005 inches at eighteen axial stations 
upstream of, within, and downstream of the interaction region. Miniature 
total head tubes, flattened to a dimension of 0.009 inches high by 0,026 inches 
wide were used for tlie pitot profiles. Velocity profiles upstream of the 
first incident shock wave, between the first reflected and second incident 
shock waves, and downstream of the second reflected shock wave were calculated 
from the pitot profiles assuming isoenergetic flow. A calibrated venturi 
meter was used to measure the bleed flow rate. 

ANALYSIS 

Figure 2 shows the flow model used in the analysis. R is the radial 
distance from the tunnel centerline and x is the distance do\-mstream from 
the tip of the generator. Conditions at the Station 1 are assumed to be 
kno\m. The object of the analysis, given the shock generator shape and 
position, is to compute the locations of the reflected shock waves and the 
boundary layer properties at successive stations along the wall. 

The computation is carried out in three steps associated with three 
subregions into which the boundary layer flow is divided. 

1) Region I extends from station 1 where the first incident shock wave 
reaches the boundary layer edge to station 3 where the reflected shock wave 
emerges from the boundary layer. Surface 2 is the stream surface which passes 
through the intersection of the incident conical shock with the boundary layer 
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edge. The location of tiis surface and the pressure distribution along it 
are obtained from an inviscld conical flow solution. Using this surface, 
the tunnel wall, and the planes normal to the wall at 1 and 3 to define 
a control volume, a control volume analysis of the region may be used to 
determine the length of the interaction, and the boundary layer thickness 
and shape at 3. In the analysis by Seebaugh, Paynter and Childs^ the method 
of solution was as follows. The velocity profile at station 1 was repre- 
sented by a least-squares-fit power-law profile. The downstream boundary 
layer thickness and power-law profile were then obtained by solving the con- 
tinuity end x-momentum equations for the control volume. The interaction 
length, x^hich is related to 6^ by the geometry of the edge stream surface, 
was determined in the solution. In the analysis, the pressure at 3 was taken 
to be the wall static pressure at station 3 as determined from an inviscid 
flow solution. The wall shear force was neglected and mass entrainment into 

the boundary layer through the interaction was ignored. 

2 3 

In a subsequent analysis Mathews ’ used a law-of-the-wall/lav7-of-the- 

wake profile to replace the power-law profile for representing the velocity 

distribution upstream and downstream of the interaction. In Mathews’ analysis 

the value of at station 1 was known from a least-squares fit of the wall- 

wake profile to the upstream velocity profile, a tentative value of C_ was 

■^3 

determined by the control volume- analysis. This made it possible to take the 

shear force on tlie wall into consideration in an approximate way, l.e. , by using 

an average of C_ and Cp to give an average wall shear stress. This scheme 
^1 ^3 

should be reasonably good for weak interactions but will overpredict the effect 
of wall shear as the shock strength is increased, and a condition of flow se- 
paration is approached. An iterative technique led to a final value for Cp 
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and an average wall shear stress. Mathews also allowed for boundary layer 
mass entrainment in the interaction region in an approximate way by taking it 
to be equal to the entrainment rate ahead of the interaction. Ilis experimental 
results supported this assimiption. Mathews* analysis led to somewhat 
better agreement betv;een predicted and measured downstream boundary properties 
than had been obtained in the Seebaugh analysis. 

In the analysis used to predict the results reported here, additional modi- 
fications to the control volume analysis have been made. Recently, Sun and 
Childs^ have developed an improved v/all-wake velocity profile for turbulent 
isoenergetic compressible boundary layer flow. The modified profile may be 
expressed in the form 


u (B^ + . f , 2A^ - B wi I ^ r 1 n 

— - A sin (arc sln( — r [1 + - -r ( In n 


. 2 avl/2 2 . M a. /I avl/2v . 

+ — (1-n)' ---ln(i+(i-n) )) 

o Si 


TI VI 

I_L n 

K u* ^ 
e 


+ cos hit)]} + 


2A^ 


( 1 ) 


or 


U 1 ^ r ♦ _l/2 n 1 1 , 2 ,, a^l/2 

— = \ sin {arc sin a [1 + - --x (In n + — (1 - n ) 

u 1/2 ^ u* a 

e a e 


— In (1 + (1 + ~ (1 + cos mr)]} (2) 

a IV u** 

e 


where 


n/K = (1/2) {(u*/u ) - (1/K) In (6u /v ) - 5.1 + 0.6M/aK} (3) 

ex T w 
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For a->«, Eq. 2 reduces to the velocity profile used by Mathews in his 
analysis. While Mathews' profile has been found to provide a good repre- 
sentation of turbulent compressible boundary layer velocity profiles, it 
does have the shortcoming that the velocity gradient does not go to zero at 
y “ With a = 1, on the other hand, this shortcoming is avoided. Furthermore, 
more realistic values of 5 are found when a “ 1. In the results reported 
here the velocity profiles have been assumed to be given the modified wall- 
wake profile with a *= 1. 

One other variation on the earlier methods has been incorporated into 
the present analysis. Namely, the flow direction in the boundary layer down- 
stream of the interaction has been taken as the average of the value at the 
wall (i.e., zero) and at y = as determined from the inviscid solution. 

The pressure at 3 has then been based on the average flow direction. In the 
analyses by Seebaugh and Mathews the flow direction downstream of an inter- 
action was taken to be parallel to the wall. o 

Figures 3a and 3b show the control volumes used in the present analysis. 
Although boundary layer bleed was not employed at the^ first interaction site, 
the control volumes shown do allow for that possibility. For the control 
volumes, shown, the continuity equation may be expressed in the form: 


2 IT p u R d R 


2npuRdR + BL 


while the x-momentum equation may be written as: 
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2 R d R - 


R -6, 

w 3 

2 IT P^ R d R - 
w 1 E 


P. 2if R d R “ 2 IT R Lx 
3 w w 

R -5. 
w 3 


[®w 2 ' T" 2 

«j 2wpu RdR- 2npu RdR+I, 


R -6^ 
w 3 


R “6,-A_ 
w 1 £ 


where x = (x + x )/2 and P„ Is the average static pressure over the 
w «3 3 

boundary layer at 3. With suitable representation of the magnitude of 

X 

and with the assumption that the velocity profile at 3 is given by Equation (2), 
the equations may be solved for L, 6^ and the profile shape at 3. Comparable 
equations may be written for the second interaction site. 

In the control volume equations given above the conditions at the beginning 
of the interaction, as well as the mass bleed rate, m„, are assumed to be known. 

ij 

Conditions along surface 2 are determined from an inviscid flow analysis for 

the given tunnel and shock generator geometry. 

The x-momentum of the bleed flow depends on the manner in which the bleed 

2 

flow is accomplished. In the analyses by Seebaugh and Mathet^s computations 

were made for three bleed models: porous-wall suctioh, slot suction and scoop 

suction* Figure 3a shows the porous wall model. With this model, the x-momen~ 

turn of the bleed flow, I- , was assumed to be zero. Figure 3b shov7s the 

a 

X 

slot-suction model. With slot suction, I_ was assumed to have the same value 

i5 

X 

as that possessed by the bleed mass as it entered the control volume, i.e. , 


2 ir p u R d R 


where R„ is determined from 


2 IT p u* R d R 
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2) Region II extends from station 3 to station 5. This is the region of 
boundary layer flow between the first reflected shock and the second incident 
shock* Since no shock interactions are present in this region, the axial pressure 
gradient is relatively low. Starting with conditions at 3 as determined by the 
control volume analysis of the first interaction, changes in the boundary layer 
properties to station 5 are computed using a turbulent boundary layer program 
suggested by Paynter and Schuehle • The program uses a wall— wake profile to 
represent the velocity profile and the entrainment function concept proposed by 
Green^ to solve the boundary layer equations. An inviscid flow solution is 
used to provide the wall static pressure distribution needed for the boundary 
layer solution in this region. The inviscid solution, however, was obtained in 
a manner which allowed for the effects of the first shock interaction. The meth- 
od employed was to use an artificial wall position in the interaction region 
which would cause the reflected inviscid shock wave position to match that deter- 
mined by the control volume analysis, 

' 3) Region III extends from station 5 to station 7. This region covers the 

the interaction of the second Incident shock wave with the boundary layer. The 

method here is similar to that for Region I, except that the flow at control 

surface 6 is no longer conical. Conditions along this surface were obtained from 

a method of characteristics solution for flow past the double-cone centerbody. 

In the characteristics solution the interaction of the first reflected shock 

\j^rh the second incident shock must be considered. The location in the flow 

field of the reflected shock was determined using the artificial wall position, 

7 

In a recent analysis of inlet flow fields by Reyhner and Hickox the effect of 
the shock wave interaction on the Inviscid flow was taken into account 
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by first obtaining a control volume solution for the boundary layer properties 
downstream of the interaction. Then, using an effective surface defined by 
the boundary layer displacement thicknesses upstream and downstream of the 
interaction and using a patching technique across the interaction region to 
construct an effective displacement surface for that region, the inviscid flow 
solution was obtained for the effective surface. A comparable technique was 
tried in the work reported here but was not as successful as the scheme of 
using the simple reflection off the artificial wall. 

For the Mach 2.82 flow with bleed, the bleed models described under Step 
1) and depicted schematically in Figures 3a and 3b were used. The bleed holes 
used in the study were drilled normal to the wind tunnel wall. At first thought, 
then, the porous *“wall model might appear to provide a better representation of 
the bleed flow. However, the bleed-hole diameter of 0.052 inches was on the 
order of one— half the boundary layer thickness at station 5. Thus, there should 
be x—iu omentum associated with the bleed flow and the bleed flow behavior might 
then be expected to lie somewhere between that for porous— wall suction and slot 
suction. As will be discussed in the section on results, this appears to have 
been the case. 


RESULTS 

The results of the analysis are shown in Figures 4 through 11. Comparisons 
are made between predicted and measured results. The data for the *=2.82 

g 

flow are from an investigation by Teeter. 

Figure 4 shows comparisons of the experimental and predicted shock wave 
patterns and boundary layer thicknesses for the Mach 3.82 flow. Also shown 
is the pressure distribution at the tunnel side wall as a function of the 
distance aft of the cone tip. The triangular points shown for the analysis 
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in the static pressure plot were determined by using the artificial wall 
and the inviscid flow solution. The predicted and observed values are seen 
to be in good agreement along the entire length of the double-shock interaction. 

Figure 5 shows d*, 6 and at several stations along the tunnel side wall. 
Here, also, predicted and observed results are in good agreement. The experi- 
mental values shown for have been determined by a least-squares fit of the 
modified wall-wake velocity profile to the experimentally determined velocity 
profiles . 

Figure 6 shows Mach Number profiles for the boundary layer at the down- 
stream end of the second interaction. The analysis predicts the end of the 
second interaction to be at x = 3.75 inches. Since profiles, were not taken 
at this specific station, profiles taken just upstream, x - 3.70 inches, and 
just downstream, x = 3-80 inches, are shown for comparison. It is apparent 
that the analysis leads to a profile which provides a good representation 
of the experimentally determined profiles near the interaction end* 

Figure 7 shows boundary layer thicknesses, shock wave patterns and wall 
static pressure distributions for the Mach 2.82 flow with 2.8 percent boundary 
layer bleed at the second interaction site. Figure 8 shows 6*, 6 and for 
this flow. ^ 

Two sets of predicted results are given, one for the porous-wall suction 
model, the other for slot suction. As is shown, tlie differences between tlie 
results for the two suction models are not large. Differences in predicted 
results with the two models are due solely to the differences in values assigned 
to the x-momentum of the bleed flux. Since the bleed rate is low, the x-momen- 
tum associated with the slot suction model is small and not too different 
from the zero value for porous suction. The predicted and measured results 
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are In reasonably good agreement. 

Figure 9 shows boundary layer thicknesses, shock wave patterns and wall 
static pressure distributions for ** 2.82 with 5.0 percent bleed. Values 
for 6*, e and ar^ shown in Figure 10, while Mach Number profiles downstream 
of the second interaction are shown in Figure 11. The flow conditions up to 
the second interaction are the same as those for the flow with 2.8 percent 
bleed. With the higher bleed rate the difference between the results for 
porous wall and slot suction are much more pronounced than with 2.8 percent 
bleed. The slot-suction model gives a reflected shock location which is in 
better agreement with the observed results. On the other hand, the values 
of 6*, 0 and obtained with the porous wall model agree better with experi- 
mental values than do the slot suction results. As was pointed out in the 
section on analysis, the Heed hole diameter of 0.052 inches was on the order 
of one-half of the boundary layer thickness so that the bleed flow behavior 
might be expected to lie between that for porous wall suction and slot suction 

The x-momeatum of the bleed value might then, in turn, be expected to lie 

\ 

between the values used with the two models. Indeed, the use of a bleed 
flow momentum flux between the two limits would lead to better overall agree- 
ment between predicted and measured values of 6*, 6 and Even then, 

however, the predicted interaction length would be too long. It should be 
remarked that the slot— suction bleed flow model does not allow for the tur- 
bulent shear stress along the stream surface separating the bleed flow from 
the main body of the flow, nor for tlie wall shear stress. Nor is the pressure 

force along the separating stream surface considered in estimating . 

X 

Tlie effects of the pressure force and wall shear tend to cancel the effect 
of shear on the separating stream surface, but the extent to which they do so 
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is not known. Further study is needed on the details of the bleed flow 
behavior, including the effect of bleed-hole bluntness, before the effects 
of bleed configuration on the shock, wave interaction can be resolved. 

CONCLUSIONS 

A. control volume analysis method, employed in conjunction with a tur- 
bulent boundary layer computation scheme, has been used to predict the flow 
field downstream of successive shock wave boundary layer interactions for 

flows at M = 3.82 and 2.82. The effects of boundary layer bleed at the 
00 

second interaction site have been considered. For flow with low bleed rates 
or no bleed the predicted interaction lengths, and wall static pressures, 
as well as the boundary layer properties downstream of the interactions show 
good agreement with measured results. As the bleed flow rate is increased, 
differences between predicted and measured results increase, indicating a 
need for improvements in the bleed flow model. 
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